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We use the hydrodynamic representation of the Gross -Pitaevskii/Nonlinear Schrodinger equation 
in order to analyze the dynamics of macroscopic tunneling process. We observe a tendency to a wave 
breaking and shock formation during the early stages of the tunneling process. A blip in the density 
distribution appears in the outskirts of the barrier and under proper conditions it may transform 
into a bright soliton. Our approach, based on the theory of shock formation in solutions of Burgers 
equation, allows us to find the parameters of the ejected blip (or soliton if formed) including the 
velocity of its propagation. The blip in the density is formed regardless of the value and sign of 
the nonlinearity parameter. However a soliton may be formed only if this parameter is negative 
(attraction) and large enough. A criterion is proposed. An ejection of a soliton is also observed 
numerically. We demonstrate, theoretically and numerically, controlled formation of soliton through 
tunneling. The mass of the ejected soliton is controlled by the initial state. 

PACS numbers: 82.20.Xr, 03.75.Kk, 05.90.+m 



I. INTRODUCTION 



The recent advances in the experiments on real Bose - Einstein condensates(BEC)i2ii^ and nonlinear optical waves^ 
have generated a huge body of works on the theoretical side based on the Gross-Pitaevskii (GV)^ (see also Refs. 
ISiSTO.) or Nonlinear Schrodinger (NLS) equation^'^^ 
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The dynamics of solutions of this equation are very complex and rich. The phenomena of coherence 
scopic tunnelingiiii^, vortex formation )^^d^d^i^^i^°'^^ instabilities, focusing and blowup are all concepts related to the 
nonlinear nature of the systems. Dark solitons, or kink-wise states, i.e. states with dynamically stable propagating 
density minima, are expected in condensates with repulsive interactions. They have been predicted for one dimen- 
sional BEG22i22ii2ii24 and may occur in higher dimensions as well. There have been several suggestions for techniques 
to engineer dark solitons in BEG^^ and these were indeed successfully created and observed^^'^^. The two groups used 
far off resonance laser beam pulse to shift the matter wave phase thus creating density minima. 

As for bright solitons, the situation is more subtle, since in such systems instability of the gas is unavoidable above 
a critical particle number where the zero point kinetic energy does not suffice to balance the collapse mechanism. The 
latter renders the system to behaviors, which deviate significantly from the mean field validity range. The critical 
number has been calculated theoretically for ^Li to be Ncr — 1400 ^^1'^'-' , which is consistent with the experimental 
measurements. A condensate with a limited number of atoms, however, can be stabilized by confinement in a trap. 
In 1995, the first evidences for a BEG in ^Li atomic gas with attractive interactions were reported'^^. Later that 
year, the stability of solitons created in condensates in a small harmonic trap, constrained to one dimensional motion, 
was predicted from numerical calculations^S. For the 3d case, however, the solitons were predicted to be stable for 
modest ranges of nonlinearities. In 1998 an analytic solution was obtained for a BEG bright soliton creation in an 
asymmetric, cigar shaped trap. It was then shown that the solitons do not expand when confinement in one direction 
is lifted-- . A major stepping stone was undoubtedly the production and observation by two experimental groups^*^ 
of bright solitons in ultra cold ''Li gas, released from a one dimensional trap. Both groups reported propagation 
without dispersion over macroscopic distances. The latter group also observed a propagation of a soliton train. 

The question of BEG bright soliton formation, stability and dynamics is far from being solved and the ongoing 
research is very active. Ways to stabilize two dimensional BEG solitons by using spatial modulation of the interac- 
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tion strength^ and by using rapid oscillations between repulsive and attractive interactions^^- were suggested. The 
existence of vortex solitons in periodic potentials (optical lattices) was revealed in Refs. ,3&.39u40l Exact solutions for 
the dynamics of 1-D trapped BEC bright solitons with a time dependent interaction strength were found in Ref. l4ll . 
Interference of tunneling BEC matter waves in an optical array was observed**^ and existence of bounce solutions in 
macroscopic tunneling was investigated in Ref. |43, . Dynamics of an initially nonuniform bright soliton was studied 
in^ as was friction and diffusion of bright solitons by the thermal cloud"'^. Dynamics of a bright BEC soliton in 
an expulsive potential were investigated in^. Finally, pulsed macroscopic quantum tunneling of BEC are expected, 
which are induced by scattering of soliton on Guassian potential barrier™. 

Remarkably, many similar phenomena are observed in light propagation, when we have to turn to electromag- 
netic wave propagation and penetration into media with different refraction indices instead of matter wave dynamics. 
Atomic correlations now correspond to the coherence of laser light, while many-body (mean- field) interactions corre- 
spond to the Kerr nonlinearity. Examples of parallel dynamics include soliton formation and modulation instability— 
in the focusing (attractive) case and dispersive shock waves in the defocusing (repulsive) ca.seM^ Multi-species 
condensates relate to multi-component, or vector, beams of light, while periodic potentials for both the atomic and 
photonic systems have been demonstrated using standing light waves.^^^^ Of course, there are also significant dif- 
ferences between the two systems, particularly when atomic excitations and quantum (vs. classical) statistics are 
involved. In these cases, too, it is useful to contrast optics with BEC in order to better understand the underlying 
dynamics of both. 

Most of the theoretical analysis of the works mentioned above have so far been dealt with by a combination of 
numerical schemes (e.g. Ref. [l3lfT3 ) and finite dimensional phenomenological models. Furthermore, all assume long 
time existence of the solitons in the gas. None has considered the problem of tracing the mechanisms responsible for 
the actual formation of the soliton in the course of tunneling. In other words, the problem of the short time dynamics 
of tunneling has so far not been addressed. The tendency to shock formation and creation of soliton at the early 
stages of the tunneling process is another new aspect of the theory presented in this paper. As we shall see below an 
interesting aspect of this process is a possibility to control such soliton parameters as mass, geometrical factors and 
velocity. 

The hydrodynamic formulation for the Schrodinger equation was originally proposed in Ref. 50. A similar approach 
is also well-known in the linear and nonlinear optics (see, e.gj ^^i^^ ). Recently hydrodynamic formalism received much 
attention^i^i^i^iS as a useful tool for analyzing the GP/NLS equation. The time independent problem has been 
studied as far back as the early 1950s. We apply this approach to one-dimensional systems. Its generalization to 
higher dimensions is straightforward although may require a special consideration of vortices. 

We study nonlinear phenomena in macroscopic tunneling of a BEC gas or optical systems. Employing the hydro- 
dynamic representation we analyze the time dependent GP /NLS equation ^ and obtain the dynamics of a trapped 
droplet tunneling through a barrier both on the short and long time scales. We predict a splitting process, in which a 
blip in the density distribution is formed at short times outside the confining potential. We find the conditions, under 
which it may evolve into an outgoing bright soliton. Our approach allows for an analytical calculation of its param- 
eters including the velocity and energy. We also show a numerical evidence for the blip and soliton formation. This 
theory allows one to design a structure, in which we can fully control the parameters of the ejected soliton, including 
its velocity and mass fraction split off of the initially trapped BEC. The latter observation also indicates a way to 
extract a stable BEC soliton out of a less stable one. These are feasible processes. Their experimental implementation 
may be carried out, e.g. by measuring light propagation in samples with spatially modulated refraction inde x^^i^" . 

II. GENERAL APPROACH 

Eq. ([T]) for a complex wave function ^'(a;,i) is equivalently written as two equations for two real functions: the 
continuity equation 
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Here 

= (4) 

is quantum potential (called also quantum pressure) (QP). The effective potential Veff{x,t) is a functional of the 
density distribution p{x, t) and may vary in time so that generally the two equation should be analyzed together. 

Sometimes we have information about the density distribution and may solve Eq. ([3]) separately. The diffusion 
term 

in the right hand side of Eq. ^ is introduced as a singular perturbation^^ which ensures its dissipative regularization, 



^-v{x,t)+v{x,t)^^v{x,t) = ^§-,v{x,t) - ^§^V.ffix,t). 



(5) 



We expect a tendency to formation of a shock wave in the velocity and density distributions. It is well know that 
the Burgers equation with the diffusion term has step-like solutions, which survive even if the 'diffusion' coefhcient C 
tends to zero. However this type of a solution does not appear if C = from the very beginning (see e.g., Eqs. (4.1) 
and (4.2) in Ref. IstI ). Therefore we first analyze Eq. ([5]) for finite values of the coefficient ^ and then take the limit 
C — ^ in the final results. 

We may now apply the Cole-Hopf transformation 



v{x, t) 



C '^x{x,t) 



m (p(x,t) 

so that the new function ip{x,t) satisfies the linear diffusion equation with a source 

V't{x,t) = -^ip^^{x,t) + ^Veffix,t)(p{x,t). 

The Green function of Eq. ^ can be represented by the Wiener path integral 
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has the form of an action for a particle with the mass m moving along a path in the potential 14/ / . Distinguishing 
as usually the contribution of the classical path Xc{t) we get the saddle point approximation for the Green function. 
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This approximation leads in fact to exact results in the limit C ^ 0. In particular the pre-exponential factor F{t) 
plays no role in this limit. We look for the solution (p{x, t) of Eq. ((T]), which at < = has the form 



ipo{x) 
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Returning to the Cole-Hopf transformation ([6]) we understand that the function St^{x) must be chosen to satisfy the 
condition 



Vq{x) 



1 dSo{x) 
m dx 



where vq{x) is the initial velocity field at t = of the original physical problem. 
The solution of Eq. ^ with this initial condition then reads 
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ip{x,t)= I dx'G{x,x',t)(po{x') = F{t) / dx'e ? (10) 



The integration in Eq. (fTU|) is carried out around the saddle point defined by 
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where x{x, t) is obtained by solving Eq. (Ilip with respect to x' for given x and t. Substituting solution (I12p into Eq. 
using the initial condition (lll|l and taking the limit C ~^ one gets the velocity field in the form 

9S'(a;,S(a;,t),f) 
= 9^ 

where the last equality follows from Eq. (|lip . 

This result has a simple interpretation. As mentioned above the function S{x,x,t) is the mechanical action of a 
particle (to be called 'tracer' below) with the mass m moving from the point a; to a; during the time t in the potential 
Veff. vo{x) is its initial velocity. In other words we have to solve the following equivalent problem. At a time t an 
observer measures the velocity field in an 'effective fluid' flowing in the external potential Vef f{x, t) under the condition 
that the initial velocity field is {x) . The possible compressibility of the fluid is accounted for by the dependence of 
the quantum pressure and T4//(x) on the density distribution. If the measurement is carried out in the point x the 
observer sees the tracer, which has started at the time from a point x with the velocity vq{x). The question is what 
the tracer's velocity measured by the observer is. Solving this problem we obtain the velocity distribution of the fluid 
flowing in the effective potential (fT5|) , which according to the above procedure coincides with the velocity field of the 
actual quantum fiuid. Below we analyze this problem for a model choice of the effective potential Veff and obtain 
the velocity field as a function of time t. 

III. MODEL 
A. Adiabatic approximation 

We consider here tunneling escape from the one-dimensional potential trap, Fig. [T] Tunneling of a particle with 



i^xt(x)> 



Vo(x) 




FIG. 1: The right half view of the trap potential Vext{x) keeping the droplet. The dashed line shows the auxiliary potential 
Vo{x) which differs from Vext{x) only to the right from the top of the potential barrier. The dashed-dotted hne corresponds to 
the chemical potential /i for a given number TV of particles. 
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the energy fi from the trap is characterized by the time of tunneling escape 



Ttun ^ V 

where v is the frequency of oscillations within the trap, and 

I rb 

A = - / dx^/2m{Vext{x) - y) 

is the tunneling integral of the classically forbidden underbarrier region {a,b). There is also another important 
characteristic of tunneling which is the time Ttr needed to traverse the undcr-barrier rcgion^^. This duality in 
tunneling characterization was first introduced in Ref. [H and has been dealt with in many works since^^'^'*'^^-^^. 
The traversal time is a matter of intensive research and attracts attention both from fundamental and technological 
reasons. Among many possible definitions of the traversal time the 'semiclassical' traversal time 

m r'' . 1 



proposed in Ref. '66' is most popular and will be used below. It is obtained by modulating the potential barrier with a 
small cosine perturbation and finding the traversal time as the crossover between high and low frequency behaviors. 

In what follows we embrace the definition (|13p for the traversal time and rely on the additional important observation 
that the density field varies on the time scale Tt„„, whereas the velocity field varies on the timescale Ttr^- The 
inequality Ttr ^ Ttum valid for typical barriers, assumes dynamics of the velocity field at virtually constant density. 
We may therefore use this assumption in order to apply the adiabatic approximation in the first iteration. The latter 
will end, by definition, when changes in the velocity field will become significant. The second step will be to calculate 
the time dependence of the density, p{x,t), for the time dependent velocity field found in the first iteration. 

We consider now a droplet trapped in the potential Vextix), from which it escapes due to tunneling through the 
potential barrier. In order to define the initial density distribution po{x) we introduce an auxiliary confining potential 
Vo{x), from which tunneling is impossible and a stationary state is formed. This potential coincides with Vext{x) for 
small X up to the top of the potential barrier but differs for larger x (dashed line in fig. [T]). The density distribution 
Pq{x] N) of the stationary state in the auxiliary potential Vo{x) is obtained from Eq. ^ at v{x, t) — 0, which becomes 



^Po{x-N) = pi{N) ^pa{x-N). (14) 



Eq. (jl4p is in fact the stationary GP/NLS equation with the potential Vo{x) and determines therefore the initial 
density distribution po{x;N). Here fJ.{N) is the chemical potential. The calculations is carried out for a given total 
number N of the particles in the trap. 

Using Eq. p4|) for po{x]N) we get that the effective potential in Eq. ([5]) becomes for a given number N{t) of 
particles becomes 



Veff{x) = Vextix) - Voix) + + V,u{{p{x, t)}) - V,u{{po{x; N (t))}) + X[p{x, t) - po{x; Nit))]. (15) 

The number of particles N{t) in the trap may vary slowly with time. Generally the effective potential follows the 
variation of density distribution. However, the adiabatic approximation implies that the density field varies much 
slower than the velocity field. It ensures also that the QP follows also the change of the total number of particles in 
the trap in the course of time, which is reflected by a slow variation of the chemical potential p, with slowly changing 
number of particles A^(t).— Hence we may, to within a good precision, assume that p{x, t) = p[j{x] N{t)) in Eq. (|15p . 
meaning that 

V,ff{x) = Vextix) - V^{x) - Uo. (16) 

where Uq is the asymptotic value of the difference Vextix) — Voix) at large x. The chemical potential may be dropped 
under the derivative over the coordinate x in Eq. ^ and we may introduce for the sake of convenience the constant 
—Uq so that Veffix) ^ at a; ^ oo. Therefore the effective potential is the shifted by Uq time independent difference 
of the actual potential forming the trap and the auxiliary potential used to form the initial density distribution. 
Then it is convenient to chose the effective potential (IT6l) in the form 



(17) 
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Since we start from the stationary state po{x; N) as the initial density distribution, we also assume the zero initial 
velocity field, Vo{x) — 0. The effective potential PT|) is time-independent so that we deal with an 'adiabatically' 
incompressible flow. The compressibility, i.e. rapid variations of p{x,t) and, hence, Vg// with time can be accounted 
for in higher iterations. 

According to the analysis in Section |TT]wc have to calculate the velocity field v{x,t) of the equivalent fluid at a given 
point X and time t with the initial conditions described above. For this sake we consider a fluid tracer which starts 
moving at a point x at the initial time t = with the initial velocity v{x, 0) = 0. The initial energy of the tracer 

e-=-^. (18) 
cosh ax 

is conserved in the Lagrange coordinates, when we follow the tracer along the path of its motion. Hence at the point 
of observation x (which is a static Euler coordinate) and at the time of observation t it becomes 

e-=eix,t)^^^ + !!^. (19) 
cosh ax ^ 

Eqs. (US]) and allow us to connect the starting point x of the tracer motion with the velocity field measured at 
the time t at the point x. 

The tracer, which has started from the point x with the energy e, reaches the point of observation x after the time 

t = I ^ dx. (20) 

(7n 




cosh^ ax 



This integral together with the initial condition (|18p implicitly defines the energy e{x,t) of the tracer observed in the 
point X at the time t. Calculating the integral in Eq. (|20p we arrive at equation 



F(w;g,r) = 2- W + 2V1 - w- wexp-^ ^ = 0. (21) 

I Vzz; sinh^ ? + 1 J 



where 



w 



e sinh^ ^ 



with < w < 1. Here we use the dimensionless time, t = ty^SUoa'^ /m, and space, ^ = ax, coordinates. It is worth 
noting that the time scale y^m/2Uoa'^ is of the order of the traversal tunneling time ()13p and appears as a natural 
scale for the time variation of the velocity field in the course of tunneling. 



B. Velocity field 

Eq. (|2T|) is solved with respect to w, so that we get the quantity e{^, r) at given ^ and r. Then we obtain the 
velocity field 

V ^ 

Eq. (I21|) is nonlinear and may have more than one solution. One can find the critical time Tc = 5.55 and position 
= 2.005 from the condition that the function F(w;^,t) in Eq. (|2ip zeros simultaneously with its first and second 

derivatives. Then for r < Tc, Eq. (|2T|) has only one solution at each value of the coordinate ^. At longer times r > Tc, 

there is a finite range of ^ values at ^ > £,c, where Eq. ((2T|) has three solutions. 

The appearance of three solutions corresponds to breakdown of the wave (see, e.g. discussion in Ref. [67t ) and to 

formation of a shock wave. In the critical region the procedure as outlined in Section [IT] should be amended. It means 

in fact that Eq. (fTTj) has several solutions and, hence, the integrand in Eq. pO)) has several saddle points. The saddle 

point, at which the action S{x,x',t) is the smallest will determine the actual velocity field in the limit C — > 0. 

The velocity field at r < Tc is shown in Fig. [21 A blip is formed on the outskirts of the barrier, which then 

approaches the wave breakdown at r = Tc Fig. [3) This type of behavior is characteristic of the shock wave formation. 
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FIG. 2: Velocity profiles (in the yj2U(,/m units), for increasing times, r — 1,2,3,4,5. 
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FIG. 3: Velocity profile (in the \/2Uo /m units) at r = Tc where breakdown of the wave occur. 



We do not follow the further evolution of the velocity field since we have stretched the adiabatic approximation to 
its limit. In fact we should consider this result as the first iteration and calculate the evolution of the density field 
by means of the continuity equation ([2]) , which may introduce corrections to the effective potential (|17|) in the region 
where the shock wave is going to be formed. 

Although the fully developed shock with a sharp step in the density distribution would be an artifact of the 
approximation of incompressible flow used to calculate the velocity field, its analysis allows for an estimate of the 
speed of blip propagation. Solving simultaneously equations 



at a given time t > we find two values of the coordinate ^ between which the shock occurs. Using the fact that w 
is small in this region, we obtain the upper limit of the blip velocity 



It is worth emphasizing that Vf, is of the order of the velocity, with which the tunneling particles traverse the classically 
forbidden barrier region. 



The development of a blip in the velocity field results in a local increase of the density at t < Tc (see Fig. U) . The 
variation of p{x, t) near the blip is found from the continuity equation ([2]) assuming similarly to Ref. [s^ that in the 
region where the blip is formed (i.e. outside the trap) the initial density distribution is 



F{w] ^, r) = 0, and 



= 




C. Density field 
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FIG. 4: The blip in the density distribution at r = 1, 3, 5 for /3 — 0.25a. The lines become dashed inside the trap where the 
calculation errors may be large. The shape of the effective potential (out of scale) is shown by the dash-dotted line. 

where (3 = 2^y2m(JJo — /i) and /5 is a constant. We solve equation 

^^=-^[po(x,tMx,i)]. (22) 
where v{x,t) is the time dependent velocity field calculated in Subsection IIII Bl The relative variation of the density 



Prelim, t) = — = dr - 

Po(0 Jo 2 



a ot, 



(23) 



is obtained by integrating Eq. (1^^ over time. 

The calculated distribution PreiiS.jT) of the relative density variation is shown in Fig. 3] at several times t < r^. 
The blip in the density propagates with the velocity Vb leaving a depleted region behind. Although the depletion is 
a real effect the approximate procedure leads to a too strong depletion inside the trap as shown by the dashed lines 
meaning that we cross the applicability limits of the procedure. 

We clearly see a tendency to the wave breakdown and to formation of a shock at times approaching the critical 
time Tc. This however doesn't happen and the wave doesn't break down and form a shock in the density field p{x), 
since such a development is inhibited by the quantum potential (|4]). The latter contains the second derivative of the 
density which may become large in the area of the blip formation even though the absolute value of Prei{£,, t) may be 
still small. The corresponding correction to the effective potential (fT51) prevents the sharp shock formation. 



1. Formation of a soliton 

At r = 3 the blip moves in a weak potential in the outskirts of the trap, and we deal with a new problem of a 
packet propagating with a the velocity vt, which may or may not transform into a soliton. The direct comparison 
shows that the shape of the blip is close to that of the soliton. The latter is formed if the blip energy Euip in its 
center of mass coordinate system is negative (e.g., Ref. [H and references therein). This condition may be fulfilled 
only if the interaction parameter A = 2mX/h? is negative (attractive interaction). Calculating E^Hp by means of the 
Gross-Pitaevskii surface energy functional we get the inequality 

|A|p>(12.95a2 + 0.123^), (24) 

i.e. for each value of Uq and a there is a lower limit for the interaction strength above which a soliton may be 
formed. Its width is about twice the width of the trap and it contains about 10% of the initial packet. For the typical 
parameters of the currently available systems, m — 7 A.M.U = 11.69 • 10~^^kg (for Li atoms), Uq = 10~^'^J, a = 10* 
m^^, and p — 10^^ m^'^ we get |A| > 2.38 • 10~^m. Finally, a soliton is formed if the interaction coefficient satisfies 
the condition |A| > 0.73 • 10~'*^J • m^, which is typically fulfilled and has been measured experimentally (see, e.g.^"*). 

It is emphasized that in the case of a weaker negative (attractive), zero or even positive (repulsive) interaction 
parameter A the soliton is not formed, however the soliton-like blip in the density will be always formed and propagate 
far away from the trap (many hundreds of barrier widths) before being dispersed. 
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FIG. 5; Contour plot of the density distribution evolving with time. The figure shows formation of a blip outside the barrier 
at short times. 



D. Numerical solution 



We may observe as soliton is formed and ejected in the course of tunneling also in the direct numerical solution 
of the GP /NLS Eq. ([1]) carried out using the program " Kitty"—. Varying the trap potential in a wide range of its 
parameters we were always able to observe formation of the blip. Here we present an example of the computation 
carried out for the GP/NLS hamiltonian 

H = -\v^ + V{x)-\\^!{x,t)\'' 



where 



with the initial wave function 



vl/(a::,0) = cosh — . 

4 4 

Formation of a blip with a negative energy outside the potential barrier at short times is shown in Fig. [5] whereas 
Fig. [5] shows its propagation at longer times. One can clearly see two parallel lines showing the core of the blip, 
which slowly oscillates approaching the soliton shape. It is worth emphasizing that the characteristics of the soliton, 
obtained numerically, are quite close to those obtained in our above analysis, including the time and location of the 
blip formation and the velocity of its propagation. 



IV. CONCLUSIONS 



Using the fluid dynamics paradigm, we have analyzed the GP/NLS equation and demonstrated for the flrst time 
the phenomenon of a shock formation in the course of tunneling short time dynamics. It results in the formation 
and ejection of a blip in the density field of the outgoing mass. When the suitable conditions are fulfilled the latter 
may later on transform into a bright soliton. We propose a technique which allows for the analysis of the short time 
dynamics of tunneling processes which provides us with a new insight to the fundamental problem of macroscopic 
tunneling. 

We believe that the conclusion of this paper can be straightforwardly verified experimentally both in BEC tunneling 
or in nonlinear optics measurements. By engineering a trap with the proper parameters we control the tunneling 
process and formation of bright solitons. The blip splits from a bigger and narrower and thus less stable trapped 
packet, therefore a general approach of ensuring stability by the BEC tunneling is presented. 

Although the model considered here assumes that a BEC droplet is bound in a trap, the blip formation in the 
course of tunneling is the property of the barrier and is observed for other initial states as well, e.g., numerics show a 
similar effect in 2d - systems. 
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FIG. 6: Contour plot of the density distribution evolving with time. The figure shows propagation of the blip and its gradual 
conversion into a soliton at long times. The lower graph uses the frame, which moves with the velocity v = 1.047 and follows 
the blip. 



The model we analyzed can be implemented experimentally both for BEC and optical soliton devices. The latter 
requires a quite simple setup of an optical fiber, with a spatially modulated refraction index, constituting the trap 
and its outings. As the light propagates, a blip should be detected near the trap outings, which escapes towards the 
sides of the fiber, and develops into a soliton. By controlling the initial state, i.e. the initial mass of the BEC droplet 
or light intensity, we may obtain a 'soliton gun' for prescribed mass and velocity of the ejected solitons. Our approach 
applies to other nonlinear soliton dynamical effects, e.g. soliton slicing™ (by a potential bump) and generally soliton 
interacting with a potential. Another important example may be the dynamics of nonlinear models of fission, (see, 
e.g^). 
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